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Abstract 

Several modeling approaches can benefit from data and dynamics transformations that allow simplification, general¬ 
ization, and/or more conceptual treatment of the studied problem. In this work, we discuss the phasor representation, 
capable of mapping sinusoidal signals into values in the complex plane, which can be associated to vectors. As such, 
phasors allow visualizating and understanding amplitude and phase relationships, such as relative phase and phase 
coherence. The potential of phasors for diverse applications is illustrated with respect to AC linear circuits, filters, 
oscillators and the Kuramoto model. 


‘...this brave o’erhanging firmament, this majestical roof 
fretted with golden fire.’ 

William Shakespeare. 


1 Introduction 

To a good extent, scientific modeling as well as technolog¬ 
ical developments often benefit from approaches capable 
of simplifying, and/or generalizing and/or obtaining more 
effective and approachable representations of aspects of 
the studied phenomena. This is frequently the case with 
transformations such as the Fourier series, in which a wide 
range of signals can be expressed in standardize way, more 
specifically in terms of a sum of sines and cosines with spe¬ 
cific amplitudes, frequencies, and phases. Other examples 
include the Laplace transform, the Taylor series, the use 
logarithm to map products into additions, and the use of 
moments to express statistical properties of signals. 

In the present text, we focus on phasor representations 
of signals. Basically, this approach allows a sinusoidal 
signal of time - with amplitude A , frequency / (or the 
associated angular velocity uj = 2nf) and phase 9 - to 
be mapped into a respective complex number, the phasor 
Ae %e , where i = y/—l (e.g. m- Recall that this complex 
number can also be associated to a vector in the plane. 
Phasor representations disregard explicit representation 
of time as well as the signal angular velocity. These two 
quantities can nevertheless be taken into account by as¬ 


suming that the mapped phasors are rotating around the 
complex plane origin with angular velocity u. The con¬ 
cept of phasor is illustrated in Figure [l] 

The advantage of the phasor approach derives mostly 
from its ability to represent in an intuitive and simple way 
the phases of the signals, as well as their interrelationships 
(e.g. relative phase, phase coherence, etc.). In addition, in 
cases where all signals have fixed frequency, it is possible 
to perform several operations between the involved signals 
in an objective and simplified geometrical way (e.g. vector 
addition). 

As a consequence of these interesting properties, pha¬ 
sors are adopted and used in a wide range of scientific 
and technological areas, including but by no means being 
limited to: physics (e.g. oscillations), dynamical system, 
electrical and electronic engineering, control theory, com¬ 
plex systems, vibrations, acoustics, audio, neurosciences, 
among many others. 




Figure 1: Plot of the signal x(t) = ut + 6 (a), and its respective 
phasor representation in the complex plane (b). 

We shall start by revising some important properties of 
the sine and cosine functions, introducing the concept of 
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phasor representation, and then illustrating its potential 
with respect to some specific applications, namely AC lin¬ 
ear circuits, linear filters, individual harmonic oscillators, 
systems of harmonic oscillators, and the Kuramoto model 
(e.g. HEIS]). The latter corresponds to a particularly in¬ 
teresting approach to modeling synchronization phenom¬ 
ena found in nature including fireflies flashing (e.eg 0 ), 
heart and neuronal operation, clocks synchronization, and 
even aspects of economical and social dynamics. 

2 The Sine and Cosine Functions 

The sine and cosine trigonometric functions are particu¬ 
larly important in science, as they model several aspects 
of the real world, being also used in several technological 
applications. Though intrinsically simple, these functions 
exhibit some interesting, more subtle, features often aris¬ 
ing from their periodical behavior. Some of these prop¬ 
erties, namely those more closely related to phasors, are 
discussed in this section. 

Both the sine and cosine functions take an angle (mea¬ 
sured in radians or degrees) as argument, which we will 
call <h(£), so that we have sm($) and cos(<F), respec¬ 
tively. Their 27r-periodicity can now be expressed as 
sin(&) = sin(& + 2 tt) and = sin(^-\-2 tt). Though 

<F can take any real value, it is often interesting to con¬ 
sider that 0 < <F < 2tt. 

One of the first interesting properties to be noticed 
about the sine and cosine functions is that they are closely 
interrelated, in the sense that any of them can be obtained 
by displacing the other along the angle argument. For in¬ 
stance, we have that 


sin(<f>) = cos(<f> — 7r/2), (1) 

cos(<F) = sin(& + 7r/2). 

It follows from this property that what we discuss and 
learn about one of these functions can be immediate and 
easily transferred for the other. In this work, we will focus 
on the sine function without loss of generality. 

Several of the applications of the sine function involve 
expressing the angle argument as a function of time, 
i.e. <£>(£). Typically, this time-dependent argument is ex¬ 
pressed in terms of a respective frequency f and a phase 
6. The period of this function is T = 1//. Observe that 
the phase 6 corresponds to the value of the sine argument 
<F(£) at the initial time, i.e. <!>(£ = 0) = 0. The angular 
velocity (or frequency ) uj = 2irf is also often considered. 
So, we can define a sinusoidal oscillating signal x(t) as the 
following function of t: 

x{t) = A sin (2tt ft + 0) = A sin (cat + 0) . (2) 


One example of this type of signal x(t) is illustrated in 
Figure [lja). 

While the frequency / determines how fast the sig¬ 
nal changes, the phase 0 can be understood as a con¬ 
sequence of displacements with respect to the reference 
sinusoidal signal sin{2i:f). The frequency / (or the an¬ 
gular frequency uj) and the phase 6 can be understood 
as representing parameters of the function in Equation [2] 
Unlike variables , such as the independent-value t or the 
dependent-value x(£), parameters remain constant during 
a given realization of the function. For instance, a lin¬ 
ear pendulum will oscillate while keeping its parameters 
/ and 6 constant. 

We have that / is commonly measured in Hertz ( Hz ), 
uj in rad/s (radians per second), and 0 in rad. Note that 
degrees (°) can also be used instead of radians. 

An interesting situation is defined when the phase is a 
function of time, i.e. 

x(t) = A sin(ujt + 0{t)). (3) 

In particular, if 6{t) is an affine function of £, i.e. 0(t) = 
vt + #cb we have 


x(t) = A sin(ujt + vt + #o) = A sin(uj e t + #o), (4) 


with uj e = [uj + v\. In other words, in this situation 
the phase 6{t) behaves in an analogous way to an angular 
frequency. As an example, consider the situation in which 
we want to obtain X 2 (t) = sin(2ujt) from a signal x\ = 
sin(ujt + 6{t)). This can be easily achieved by making 

v — UJ. 

Recall that we have referred to the argument of the sine 
function as 4>(t). Now, the respective time derivative of 
this argument, given as 


d$(t) 

dt 


m 


(5) 


has unit rad/s and can be effectively understood 
as the instantaneous (angular) frequency of the signal 
Asin($(t)). For instance, if 6{t) = vt + 0q we have, as 
expected, that 4 >( t ) = d{ujt + vt + 6q)/ dt = uj + v = w e . 


3 Phasors 

It is often interesting and useful to visualize signals such 
x(t) in Equation [l8ja). We have already illustrated 
the Cartesian representation of the sine function in Fig¬ 
ure [lja). However, this type of representation can be 
highly redundant as a consequence of the 2tt periodicity 
of the sine. A more effective approach can be obtained by 
using the Euler equation to map the sine function in the 
complex plane (e.g. HaiaiE]), i.e. 
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x(t) = Asin($(t)) = Aim , (6) 

which can be associated to the time-varying complex 
variable z(t) = Ae % ^\ 

Figure [ljb) shows the phasor representation of the sig¬ 
nal x(t) = Asin( 2 Tift + 6 ) in Figure[lja) into the complex 
plane. Observe that the sine values are mapped into a 
circle of radius A, while its argument 0 < <F(£) < 2 tt cor¬ 
responds to the angle between each of the function points 
and the Real axis. It is also possible to understand the 
respective signal x(t) as being obtained by rotating the 
phasors counterclockwise with constant angular velocity 
uj while taking the respective projections of the phasor 
onto the imaginary axis. 

For the very same reason that the phasor representation 
allows the elimination of the redundancy derived from 
the sine periodicity, it also implies that we are no longer 
able to geometrically discriminate between values of $>(t) 
(i.e. angles) that are larger than 2 i r or smalller than 0. 
In other words, we can only visualize the function x(t) 
as defined for 0 <<£>(£)< 2 i r. The discussed Cartesian 
and complex representations of x(t) can be understood 
as being unwrapped and wrapped around the period 0 < 
<F(£) < 27r, respectively. 

Now, the smallest difference between the two phases, 
which provides a measurement of the similarity between 
the two respective phases, can be immediately appreci¬ 
ated in terms of the smallest angle between the two pha¬ 
sors in the respective complex-plane representation. Note 
that, in the general case, this angle cannot be always ob¬ 
tained mathematically by making S = O 2 — # 1 , as a conse¬ 
quence of the periodicity of the phasor representation. For 
instance, in case 0\ = 5° and 62 = 355°, we would have 
S = 350° instead of the expected smallesr angle differ¬ 
ence of just 10°. For such reasons, it becomes interesting 
to use other ways of obtaining measurements expressing 
angle differences while disregarding the possible discon¬ 
tinuities implied by the angle periodicity. An interesting 
possibility is to adopt AO = | sin (62 — 0i)|. 

4 AC Linear Circuits 

A powerful application of phasors regards the time- 
varying analysis of linear electric circuits, which is usually 
called alternated current (AC) analysis (e.g. [lj 9j), as a 
counterpoint to direct current (DC) analysis. The AC 
regime involves a transient period followed by an steady- 
state stage. Here we focus attention on the latter situa¬ 
tion, characterized by the almost complete elimination of 
the effects of the initial charges and fields in the capaci¬ 
tors and inductors. The complete behavior of the circuit 
can be obtained by linear superimposing the transient and 


Table 1: The reactances and impedances of the passive components 
R , C and L. Observe that i = e F 71- / 2 ) and —i = ed -7r / 2 ). 


component 

reactance 

impedance 

R 

R 

R 

C 

X C = ~^c 

z c = iX c = ~zjc 

L 

X L = luL 

Zl — iXL = iuL 


steady-state solutions. 

Interestingly, the linear circuit framework is not only 
very useful in electricity and electronics, but it can also 
be related and applied to other concepts and areas such as 
acoustics, vibrations, temperature, flow, hydrodynamics, 
magnetism, among others. 

AC linear circuits involve combinations of some of the 
following possible components: voltage sources (v(t), in 
Volts), current sources (i(t), in amperes, A), resistors 
(R, in ohms, fi), capacitors (C, in farads, F), and induc¬ 
tors (L, in henries, H). Sources are understood as active 
components, while resistors, capacitors and inductors are 
said to be passive. Note that, in addition to at least one 
time-varying voltage or current source, AC analysis may 
involve also other time-varying or time-constant voltage 
and/or current sources. Circuits incorporating more than 
one source can be analyzed through linear superimposition 
of the solutions found considering each of them separately. 

The time-varying signals are assumed to be in the har¬ 
monic form. For instance, we can have a voltage source 
v(t) = Asin(ojt + 0), where A E 5ft is the signal amplitude. 
This signal can be represented in the complex plane as 

v(t) = Im = AIm(e iut .e ie ). (7) 

So, the original voltage signal v(t) is transformed into 
the real number A. e luJt . e 10 . Because a linear circuit with 
a source operating at uo will have currents and voltages at 
all passive components oscillating with that same angular 
frequency, we can take into account only the amplitude 
and phase of these signals, transforming v(t) into the re¬ 
spective phasor V = Ae lG . 

In a purely resistive resistor R, the voltage and current 
across it can be related by Ohm’s law v(t) = Ri(t). This 
law can be immediately extended to capacitor and induc¬ 
tors by defining the respective reactances and impedances , 
as given in Table [l] As with resistance, the reactances and 
impedances are measure in Ohms. 

Observe that reactances assume real values, while 
impedances are purely imaginary. Let’s consider a com¬ 
ponent consisting of a resistor R in series with a capacitor 
C and an inductor L. The total reactance is defined as 

X = X C + X L (8) 
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and the total impedance is given as 

Z = R + iX = R + i(Xc + X L ) = R + T (9) 

This composed impedance device represents a general 
case of impedance in a linear AC circuit. Interestingly, 
impedances Z can be associated in a similar way to finding 
the equivalent resistance of a network of resistors. 

The above definitions allows the Ohm’s law to be 
rewritten as: 

V = ZJ, (10) 

in which / is oriented in opposition to V along the pas¬ 
sive component Z. 

Let’s consider a single capacitor C connected to a volt¬ 
age source V = Ae l ° with a generic phase 0 and ar¬ 
bitrary amplitude A. For this case, we have that I = 
V/Z = VjZ c = VK-i/uC) = iV ujC = e hV 2 ) e ^ AuC = 
l jjACe 2 ( 7r / 2+0 ) ? and consequently the voltage V lags (in the 
sense of being delayed in time) the current / by a relative 
phase 7r/2, irrespectively of a;, A , and C. The respective 
phasors are illustrated in Figure [2|a). 

Now, if instead of a capacitor, we attach an inductor 
L to the source V, we have I = V/Zl = V/(iuL) = 
-iV/(uL) = Ae^-^e 0 /(ujL) = A/{uL) e^-"/ 2 ), and 
now it is the current that lags the voltage by a relative 
phase 7r/2. The respective phasor diagram is given in 
Figur^b). 

As a more general example, consider the circuit con¬ 
taining a resistor R = 100011 and a capacitor C = 0.01 fiF 
in series with a voltage source v(t) = 2sin(2i\ft + 7r/4), 
i.e. V = 2e*( 7r / 4 ) Volts , as shown in Figure^ We adopt 
/ = 1/(2ttRC) ~ 15915.49iLz, so that c j = 10 5 rad/s , 

It follows that Xc = —l/(10 _3 )fl = —10 3 fl and Zc = 
-i 10 3 fl, so that Z = 1000-il0 3 H = 1000+10 3 e^- 7r / 2 ^. 
Now, we can determine the current I by using the complex 
Ohm’s law as: 



Figure 2: (a): In a capacitor , the phasor notation helps to visualize 
that the voltage V lags the current / by 7r/2. (b): In an inductor , 
the current / lags the voltage V by 7r/2. (c): The voltage phasors 
with respect to theexample in the text, consisting of voltage source 
V in series with a resistance R (voltage Vc) and a capacitor C 
(voltage Vl). The current phasor /, not shown, is parallel to Vr. 

5 Linear Filters 

AC circuits such as that in the RC example can be under¬ 
stood as filters (e.g. m), in the sense that they will at¬ 
tenuate the input signal (voltage or current) as a function 
of the signal frequencies uj. Filters are very important 
in several applications ranging from electronics to neu¬ 
rosciences. In this section we illustrate how the phasor 
approach can be applied to better represent and under¬ 
stand filter operation. Generally speaking, linear filters 
involve the passive components resistors, capacitors and 
inductors. A filter is often studied by considering its op¬ 
eration, i.e. how it changes the magnitude and phase of 
the input signal V in terms of the respective frequency /, 
which can be expressed as: 


I = 


V 

z 


2gh 7r / 4 ) 

1000 + 10 3 e*( -7r / 2 ) 


v / 2e i ( 7r / 2 ) =iV2mA. 


( 11 ) 

By considering R and C in terms of their separated 
impedances Zr = R = lOOOf^ and Zc = —i/{ujC) = 
—il0 3 S2, we can now obtain the voltages across these com¬ 
ponents as Vr = RI = \f2e l(<K l 2 ^ volts and Vc = ZqI = 
10 3 a/2 10 -3 = y/2 volts. Observe that we could also have 
calculated Vc from V and Vr by using Kirchhoff’s voltage 
law, i.e. V = VcAVr =>V c = V-V r = 2e i < 7r / 4 )-i V2 = 
\f2 volts. 

Note that Kirchhoff’s voltage law can be suitably rep¬ 
resented as the complex (or vector) sum in Figur^2|c). 
Observe also that the phasor / defined by the current i(t) 
is parallel to the phasor associated to Vr. 


V m 




A , pidout 

^out V 

A in Z idin 


a v e^ 0out din) = OLy e i(t)v , 


( 12 ) 


where ay is the voltage gain and fiy is the relative 
phase implied by the filter. Observe that both these gains 
are functions of the frequency /. 

An example of a linear filter corresponds to the series 
configuration of V, R and C addressed in section @ By 
allowing V to vary in terms of /, instead of fixating it 
at / = 1/(27 tRC), we have that the impedance of the 
capacitor will also change with / as Zc = — i/(ujC), so 
that I = V/Z = V/(R - i/(uC)), implying V c = Z C I = 
V/(l + iujRC ), so that 


Vc = 1 

V 1 + iujRC 


(13) 
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Figure [3] shows the phasors V, Vr and Vc as the 
frequency / increases, in steps, from /i = 500 Hz to 
fg = 50 kHz along the sense indicated by the respec¬ 
tive integer indices. Observe that, for any frequency, we 
can obtain the phasor associated to the voltage source as 
V = Vr + Vc, as illustrated by the dashed black arrows 
in this figure for fg. It can be readily verified from the 
obtained phasor representation that, for small values of 
f, Vr tends to align with V; while as / increases, Vc 
will be substantially reduced while aligning with V as a 
consequence of the reduction of the capacitive reactance 
verified for higher frequencies. 



6 Harmonic Oscillators 

Let’s now briefly consider an important physical phe¬ 
nomenon, namely harmonic oscillations. Consider the 
pendulum (e.g. [115]) in Figure [4ja). Here, a sphere (or 
bob) with mass m is attached to a pivot through a weight¬ 
less rod of length L which can rotate freely around its 
pivot. For simplicity’s sake, we will assume that all the 
mass m is concentrated at the center of mass of the sphere. 
Also, it is assumed that the pendulum is in the vacuum 
and that the pivot is frictionless, so that there is no en¬ 
ergy dissipation. The angular position /3(t) of the sphere 
is taken with respect to the vertical, and the initial an¬ 
gular position is taken as /3(0) = /3q. The initial angular 
velocity is assumed as /3(0) = 0, i.e. the sphere is initially 
at rest. 



Figure 4: Modeling the angular movement of a pendulum. A sphere 
with mass m is attached to a pivot through a weightless rod. The 
angular position (3{t) of the sphere is taken with respect to the 
horizontal axis. 


Figure 3: The behavior of the RC filter example from Section[4](R = 
100014, C = 0.01 fiF and V = 2e Z7r / 4 can be intuitively represented 
in phasor notation. 


As the magnitude of Vc attenuates with /, it becomes 
interesting to define some critical point attenuation point. 
The standard reference value is the cutoff frequency, f c , 
corresponding to where the signal power is halved. This 
implies an attenuation by ay = \[2)2 ~ 0.707 of the 
signal. It can be found that 


2t tRC' 


(14) 


For f = f c , we also have that the signal corresponding 
to V undergoes an attenuation of 20Log(\/2/2) « —3 dB 
and the output signal lags V by — 7t/4 (as is the case with 
the previous example, see Figure [2]). The passing band 
of this RC filter can the be understood as the frequency 
interval 0 < f < f c , implying this circuit to operate as a 
lowpass filter. 


The mass is under the force of gravity / = mg, which 
can be decomposed into two orthogonal components, one 
along the rod orientation, f n , and another tangential to 
the sphere displacement, f t . The action of f n is counter¬ 
acted by a respective reaction f r , as there is no possibility 
of the rod to expand or contract, so that f n + f r = 0 and 
then f r = —f n . On the other hand, the tangential com¬ 
ponent induces an angular displacement of the sphere. 
By trigonometric decomposition of /, we derive the mag¬ 
nitudes f t = f sin(/3(t)) and f n = f cos(/3(t)). Recall 
that the angular acceleration j3(t) is related to the linear 
counterpart a(t) as a(t) = Lf3(t). By Applying Newton’s 
second law, we derive: 


ft = mg sin(/3(t)) = —m aft) « — m L f3(t) =4> (15) 

P(t) = -j-sin(/3(t)). 

Though this differential equation completely specifies 
the angular displacement /3(t), it is difficult to solve as 
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a consequence of its nonlinearity. For small angular dis¬ 
placements /3(t) < 1 rad, we can adopt the approxima¬ 
tion sin(/3) ^ implying the simplified linear differential 
model 

p(t) = ‘-j,0(t)- ( l6 ) 

We observe that the name harmonic motion tradition¬ 
ally refers to the fact that the restoring force is directly 
proportional to the opposite displacement, as in the pre¬ 
vious equation. 

The solution of Equation [l6] yields /3{t) = 

/?o cos ^y/g/Lt^ = /3q sin (^^/gjLt-\- tt/2^ . Thus, 

we find that the considered pendulum undergoes a sinu¬ 
soidal oscillation around the vertical axis with angular 
amplitude of the oscillation equal to /?o- By considering 
that yJg/L = cd(t), we have / = 1/(2 n)y/g/L. By 
making 6 = tt/2 , the solution can now be expressed as 

f3(t) = /3 0 sin(ujt + 0) = /3 0 sin($(f)), (17) 

which is in the same form as Equation [2] 

Therefore, we obtain that a linear pendulum produces 
a signal that can be described by a phasor. We also have 
that <F(£) corresponds to the angular displacement of the 
phasor associated to the oscillator, i.e. e 

7 Independent Harmonic Oscilla¬ 
tors 

Let’s now consider not one, but N oscillators that are 
completely disconnected, and therefore, independent one 
another. Each of these oscillators j will have respective 
angular frequency c Oj and phase Oj. Because the phase 
is 27r-periodical, it is convenient to visualize the phase in 
the complex plane. Considering the Euler equation, we 
have: 

sin (ujt + Oj) = Im (e^^)) = Im (e i(w ^. e*') . 

(18) 

In case all the angular frequencies are equal , meaning 
that ujj = uj for any oscillator j, the phase difference be¬ 
tween any pair of oscillating signals Xi(t) is time-invariant, 
so we can conveniently visualize the oscillations only in 
terms of the respective quantities (i.e. a phasor) de¬ 
fined by the respective phases Oj. 

Figure [5] shows, in the complex plane, the distribution 
of N = 100 phases following a uniform distribution in the 
interval [—7r,7r). In order not to clutter the visualization 
of the phasors, each of them is represented only in terms 
of small circles of different collors. 



Figure 5: Phasor representation of the phases of N = 100 oscillators 
of the type Xj{t) = sin(ut + Oj), where the phases 0i were drawn 
with uniformly random probability in the interval [—7r,7r). The 
black dot near the center of the complex plane corresponds to the 
parameter n in Equation |19| Note that n —>• 0 as N —>• 00 . 

This type of phase distribution presents an intrinsic ‘or¬ 
der’ that can be evaluated in terms of the parameter: 

1 N 

II = re* = - (19) 

3 = 1 

where r gives the phase coherence and ip is the average 
phase. Observe also that II corresponds to the phasors 
center of mass. For uniformly random distributions such 
as that in Figure [5| it can be shown that r —>> 0 as TV —>> 00 . 

8 The Kuramoto Model 

The Kuramoto model was developed by Yoshiki Ku¬ 
ramoto along the 70’s as a simplified version of Winfree’s 
approach [2] to the synchronization of several coupled os¬ 
cillators. Kuramoto was then able to solve this non-linear 
dynamical problem in analytical fashion. 

In the Kuramoto model (e.g. @si), the dynamics per¬ 
formed by each oscillators is assumed to be governed by 
the following equation 

N 

&i(t) = Kijsin (&j(t) — &i(t)) = LJi+0i(t). ( 20 ) 

2=1 

Here, IT is a coupling matrix and is the angular 

displacement of oscillator i, with i = 1,..., N. Each os¬ 
cillator has its own natural angular velocity cd*. Observe 
that when K = 0, i.e. the oscillators are uncoupled, we 
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have = LJi, and each oscillator i will acquire an in¬ 
dependent angular velocity uji. Also important to keep 
in mind that the individual oscillators in the Kuramoto 
model are non-harmonic and non-linear. 

For each oscillator, we can write 


harmonic oscillations of a linear pendulum, a system of 
independent oscillators, and the Kuramoto model. Sev¬ 
eral other applications exist, and others will continue to 
be devised. 


4q(t) — uji + Oi(t) =>* — Wit + 0{t) + #o? (21) 


for i = 1,..., AT, and where Oq corresponds to a refer¬ 
ence phase. This formulation allows us to associate each 
oscillator to a respective phasor e ^ Wit+di ^\ while dis¬ 
regarding the reference phase 0$. In this way, we can 
summarize the dynamics performed by each oscillator to 
the function sin($i(t)). 

Let’s take into account the situation arising when each 
oscillator has a slightly different frequency cMore 
specifically, it is assumed that these values are drawn from 
a normal distribution with average (i and relatively small 
standard deviation a. For simplicity’s sake, but without 
loss of generality, we make fi = 0. Let’s also assume that 
all coupling weights are given as LQj = K/N , and that 
K is relatively small, e.g. K = 1. As all oscillators in 
the Kuramoto model have the same unitary amplitude, 
all phasors will lie along the circle of radius 1 centered at 
the origin of the complex plane. 

Figure [6] illustrates the unfolding of the angular dis¬ 
placements (calculated numerically) represented in 

terms of the respective phasors. Here, K « 1, fi = 0 and 


a = 0.1 rad. 

Recall that the angular velocities Wi of each oscillator 
are slightly different one another, remaining fixed along 
time. So, how can these oscillators synchronize in fre¬ 
quency even so? This takes place because the phase vari¬ 
ations 6i(t) implied by the governing Equation 
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effec¬ 
tively act so as to achieve similar effective angular veloc¬ 
ities Wi + 0i(t). 

An important aspect of the Kuramoto dynamics re¬ 
gards its specific stability. Depending on r and if, the syn¬ 
chronization can be stable or unstable , defining respective 
bifurcation diagrams (e.g. [311]). The Kuramoto model 
has motivated several extensions and variations, including 
the consideration of complex interconnections (e.g. m) 
on the respective dynamics (e-g. 0 ). 


9 Concluding Remarks 

Despite their relative simplicity, phasors allows us to vi¬ 
sualize and identify relationships between phases in sev¬ 
eral types of problems, including circuits, oscillators and 
complex systems. In this text, we presented the phasor 
notation at an introductory level and illustrated its po¬ 
tential for applications with respect to an RC circuit, the 
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Costa’s Didactic Texts - CDTs 


This is a Costa’s Didactic Text (CDT). CDTs 
intend to be a halfway point between a formal sci¬ 
entific article and a dissemination text in the sense 
that they: (i) explain and illustrate concepts in a 
more informal, graphical and accessible way than 
the typical scientific article; and, at the same time, 
(ii) provide more in-depth mathematical develop¬ 
ments than a more traditional dissemination work. 

It is hoped that CDTs can also provide integration 
and new insights and analogies concerning the 
reported concepts and methods. We hope these 
characteristics will contribute to making CDTs 
interesting both to beginners as well as to more 
senior researchers. 

Though CDTs are intended primarily for those 
who have some preliminary experience in the 
covered concepts, they can also be useful as 
summary of main topics and concepts to be learnt 
by other readers interested in the respective CDT 
theme. 

Each CDT focuses on a few interrelated concepts. 
Though attempting to be relatively self-contained, 
CDTs also aim at being relatively short. Links 
to related material are provided in order to 
complement the covered subjects. 

The complete set of CDTs can be found 
at: https://www.researchgate.net/project/ 

Costas-Didactic-Texts-CDTs. 
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Figure 6: Time unfolding of phasors obtained from the phases of N oscillators in the Kuramoto model assuming K = 1, /i = 0 and a = 0.1 
rad. The angular velocities were drawn from a normal distribution with null mean and unit variance, and are remain fixed along time. The 
plot titles indicate the respective time interaction, while the red X identifies the phasors center of mass, leaving the trajectory in black. 
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